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NOMENCLATURE 


bmep  Brake  mean  effective  pressure  (kPa) 

bsfc  Brake  specific  fuel  consumption  (g/kW-hr) 

CTM  Characteristic  time  model 

Da  Damkbhler  number 

dnoz  Nozzle  hole  diameter  (m) 

E  Activation  energy  (cal/gmol) 

EGR  Exhaust  gas  recirculation 

Kc  Equilibrium  constant  based  on  concentration 

k/f  Forward  rate  coefficient  of  the  ith  reaction 

k,r  Reverse  rate  coefficient  of  the  /th  reaction 

limp  Fuel  jet  impingement  distance  (m) 

Mno2  Molecular  weight  of  NO2 

m  Model  constant  incorporating  various  conversion  constants 

mf  Mass  of  fuel  per  cycle  (g) 

msf  Model  constant  for  soot  formation  process 

mso  Model  constant  for  soot  oxidation  process 

m'  '  Scale  factor 

N  Engine  speed  (rpm) 

NOxEI  Oxides  of  nitrogen  emissions  index  (g  NOx  as  N02/kg  fuel) 

P2  Cylinder  pressure  at  start  of  combustion  (atm) 

P3  Maximum  cylinder  pressure  (atm) 

PSR  Perfectly  stirred  reactor 

R  Universal  gas  constant  (1.986  cal/gmol  K  or  82.06  atm  L/kmol  K); 

correlation  coefficient 
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Re 

Ri 

Tc 

SootEI 

T2 

Tab 

T50% 

T(|i=i 

T,|,=3 

Vab 

We 

AP 

Eno 

<!> 

Pa 

Pd 

^fi,no 

Xid 

Ximp 

Xno 

Xno,3b 

Xsf 

Xsl,3b 

Xsl,sf 


Reynolds  number 

/th  reaction  in  a  mechanism 

Compression  ratio 

Soot  emissions  index  (g  soot/kg  fuel) 

Cylinder  temperature  at  start  of  combustion  (K) 

Burned  gas  temperature  at  the  end  of  combustion  (K) 

50%  fuel  recovery  temperature  (K) 

Stoichiometric,  adiabatic  flame  temperature  (K) 

Temperature  at  a  equivalence  ratio  of  3  (K) 

Cylinder  volume  at  the  end  of  combustion  (cm^) 

Weber  number 

Pressure  drop  across  fuel  injection  nozzle  (MPa) 

Fractional  extent  of  reaction  for  NO  formed  per  eddy  in  a  given  zone 
Fuel-air  equivalence  ratio 
Density  of  air 
Power  density 

Eddy  lifetime  in  NO  formation  zone 
Characteristic  time  for  ignition  delay 

Characteristic  time  for  fuel  impingement  on  cylinder  bowl  wall  (based  on 
the  penetration  correlation  of  Arai  and  Hiroyasu  (1990) 

Characteristic  kinetic  time  for  NO  formation  (s  or  ms) 

Characteristic  kinetic  time  for  NO  decomposition  (s  or  ms) 

Characteristic  kinetic  time  for  soot  formation  (s  or  ms) 

Eddy  lifetime  in  NO  decomposition  zone 

Eddy  lifetime  in  soot  formation  zone 

iii 


^sI,so 

'Cso 


Eddy  lifetime  in  soot  oxidation  zone 
Characteristic  kinetic  time  for  soot  oxidation  (s  or  ms) 
Kinematic  viscosity  of  fuel 


Subscripts: 


eq 

Equilibrium 

f 

Forward  reaction 

no 

Nitric  oxide 

ov 

Overall  equivalence  ratio 

r 

Reverse  reaction 

sf 

Soot  formation 

so 

Soot  oxidation 

(l)=l 

Stoichiometric 

00 

High-pressure  limit 

0 

Low-pressure  limit 

1 

Zone  1  conditions 

r 

Zone  r  conditions 

2  ' 

Zone  2  conditions 

3b 

End  of  combustion 
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1-1.  The  two-zone  model  for  NO  formation  and  decomposition  in  a  DI  5 
Diesel  combustion  chamber  as  postulated  by  Mellor  et  al.  (1998).  Zone 
r  is  intermediate  and  based  on  the  zone  1  pressure  and  flame 
temperature  but  the  zone  2  stoichiometry,  (|)ov,  for  the  engine. 

1-2.  The  contribution  of  the  Zeldovich  and  N2O  mechanisms  to  the  NO  9 
formation  rate  as  postulated  by  Mellor  et  al.  (1998).  The  information  in 
this  figure  is  presented  in  Arrhenius  format  for  typical  engine 
conditions.  Logarithms  are  presented  as  functions  of  inverse 
stoichiometric  flame  temperature  in  the  approximate  range  from  2200  to 
2900  K,  so  that  midrange  values  correspond  to  stoichiometric  values, 
with  pressure  as  a  parameter. 

1-3.  Engine-out  NOx  emissions  index  versus  model  parameter  (Eq.  (14))  for  17 
wall  impingement  datum  points  with  SOC  =  0.  The  legend  has  the 
format:  engine  speed  (RPM)/load  (kPa)/injection  pressure  (MPa).  Each 
set  of  data  is  obtained  via  an  EGR  sweep. 

1-4.  Engine-out  NOx  emissions  index  versus  model  parameter  (Eq.  (14))  for  18 
free  jet  datum  points  with  SOC  =  0.  The  legend  has  the  format;  engine 
speed  (RPM)/load  (kPa)/injection  pressure  (MPa).  Each  set  of  data  is 
obtained  via  an  EGR  sweep. 

1-5.  Schematic  of  a  DI  Diesel  spray  plume  with  regions  of  soot  formation  20 
and  oxidation  labeled. 

1-6.  Performance  map  for  a  naturally  aspirated  DI  Diesel  engine  (Heywood,  24 
1988). 
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1.0  CHARACTERISTIC  TIME  MODEL 


1.1  STATEMENT  OF  THE  PROBLEM  STUDIED 

The  goal  of  this  work  has  been  to  develop  characteristic  time  models  (CTMs)  for 
the  prediction  of  engine-out  NOx  and  soot  emissions  from  direct  injection  (DI)  Diesel 
engines.  The  CTM  approach  to  emission  prediction  involves  identifying  the 
characteristic  times  of  the  dominant  chemical  and  fluid  mechanic  subprocesses  in  the 
combustion  phenomenon  and,  based  on  these  times,  develop  a  correlation  between  the 
appropriate  model  and  pollutant.  From  a  computational  standpoint,  this  type  of  semi- 
empirical  correlation  is  relatively  simple.  Consequently,  the  model  is  useful  as  a 
preliminary  design  tool,  which  may  be  used  by  Diesel  engineers  to  reduce  the  cycle  time 
required  to  design,  develop  and  test  new  engines.  The  computational  simplicity  of  the 
model  also  makes  it  amenable  to  incorporation  into  a  CFD  code  or  a  phenomenological 
cycle  simulation  code  to  improve  the  emissions  prediction.  Still  another  perceived  end 
use  is  for  real-time  engine  emissions  control. 

In  previous  work  done  to  model  NOx  emissions  in  Diesel  engines,  Ahmad  and 
Plee  (1983)  showed  that  NO  kinetics  are  described  by  a  single  global  activation 
temperature.  The  global  activation  temperature  suggested  by  these  investigators  for  DI 
Diesels  was  34,300  K,  significantly  less  than  the  global  activation  temperature  that  is 
typically  associated  with  NO  formation  through  the  extended  Zeldovich  kinetics  in 
conventional  gas  turbine  combustors  (68,479  K,  Sawyer  et  al.,  1973;  67,976  K,  Tuttle  et 
al.,  1977).  Ahmad  and  Plee  (1983)  attribute  this  difference  to  NO  decomposition  within 
the  engine.  Thus,  the  model  of  Ahmad  and  Plee  (1983)  implies  that  the  importance  of 
NO  decomposition  relative  to  formation  in  Diesel  combustion  is  approximately  constant. 
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More  recent  analysis  by  Duffy  et  al.  (1996),  however,  suggests  that  a  single  global 
activation  temperature  does  not  model  NOx  emissions  from  modem  DI  Diesel  engines. 
Consequently,  the  one-step  global  model  of  Ahmad  and  Plee  (1983),  though  adequate  for 
the  engines  analyzed  in  that  investigation,  is  not  generally  applicable  to  DI  Diesel 
engines. 

The  work  of  Plee  and  coworkers  (Plee  et  al.,  1981;  Ahmad  and  Plee,  1983)  also 
used  a  single  global  activation  temperature  to  describe  the  kinetics  of  particulate 
emissions.  The  obtained  activation  temperature  represented  the  contributions  of  both 
soot  oxidation  and  formation,  and  varied  from  engine  to  engine.  The  variations  of  the 
activation  temperature  was  attributed  to  differences  in  mixing  and  thus  combustion 
chamber  geometry.  The  increased  mixing  causes  the  formation  of  soot  to  be  reduced 
along  with  better  recirculation  of  soot  particles  through  the  reaction  zone,  enhancing  soot 
oxidation  (Ahmad  and  Plee,  1983).  Although  both  chemical  processes  involved  in  soot 
emissions  are  important,  it  was  concluded  that  the  oxidation  process  is  dominant  and 
affected  more  by  flame  temperature  changes. 

In  this  work,  the  one-step  global  approach  to  modeling  NO  kinetics  in  Diesels  is 
replaced  by  a  skeletal  mechanism,  which  contains  the  NO  chemistry  relevant  to  high 
pressure  and  temperature  combustion  typical  of  DI  Diesels.  Additionally,  based  on  the 
work  done  under  this  program  as  well  as  laser-sheet  imaging  results  from  Dec  (1997),  a 
zonal  model  of  the  physics  of  the  NO  formation  and  decomposition  processes  in  Diesels 
is  postulated.  From  these  kinetic  and  physical  models  for  NO  emissions,  a  characteristic 
time  model  for  NOx  emissions  from  Diesel  engines  is  derived.  From  kinetic  times  for 
NO  formation  and  decomposition,  derived  here  from  first  principles,  the  factors  affecting 
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the  importance  of  NO  decomposition  are  evaluated.  Preliminary  empirical  analysis  to 
define  a  fluid  mechanic  time  for  NO  formation  has  also  been  performed  under  this 
program. 

For  the  soot  emissions  model,  the  global  activation  temperature  approach  still 
applies,  but  the  contribution  of  both  soot  formation  and  oxidation  will  be  represented  by 
individual  chemical  kinetic  and  mixing  times.  Also  based  on  the  work  of  Dec  (1997),  a 
two  zone  model  is  postulated  for  soot  with  the  importance  of  each  chemical  process 
dependent  on  the  spatial  location. 

1.2  SKELETAL  MECHANISM  FOR  NO  CHEMISTRY 

NOx  chemistry  can  proceed  through  the  extended  Zeldovich  mechanism  (Lavoie 
et  al.,  1970),  the  prompt  mechanism  (Fenimore,  1971),  the  N2O  mechanism  (Malte  and 
Pratt,  1974),  and  NO  formation  fi'om  fuel-bound  nitrogen.  The  last  pathway  is  negligible 
for  modem  Diesel  fuels,  which  contain  low  amounts  of  organic  nitrogen. 

To  formulate  a  skeletal  mechanism  for  NOx  emissions  in  Diesel  engines,  it  is 
necessary  to  consider  the  different  pathways  of  NO  formation  and  decomposition  and, 
through  kinetic  analysis,  determine  those  pathways  that  are  significant  at  the  operating 
conditions  of  interest.  Mellor  et  al.  (1998)  present  a  review  of  such  kinetic  analyses. 
From  this  review  it  is  concluded  that,  at  conditions  typical  of  Diesel  combustion  (i.e., 
stoichiometric  and  lean  flame  temperatures  and  high  pressures),  the  significant  paths  for 
NO  chemistry  are  the  extended  Zeldovich  (Lavoie  et  al,  1970)  and  nitrous  oxide  (Malte 
and  Pratt,  1974)  mechanisms. 

Based  on  the  review  discussed  above,  Mellor  et  al.  (1998)  postulate  a  skeletal 
mechanism  for  NO  chemistry  in  DI  Diesel  engines.  This  mechanism  is  expressed  as 
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O  +  N20NO  +  N  (Rl) 

N  +  O2ONO  +  O  (R2) 

N  +  OHoNO  +  H  (R3) 

N2O  +  O  o  2N0  (R4) 

O  +  N2  +  M0N2O  +  M  (R5) 

O2  +  M02O  +  M  (R6) 


N20  +  MoN0  +  NH...=>2N0  (R7) 

where  reactions  (Rl),  (R2),  (R3),  and  (R6)  are  from  the  extended  Zeldovich  mechanism 
and  (R4),  (R5),  and  (R7)  are  from  the  N2O  mechanism. 

Along  with  the  kinetic  studies  reviewed  by  Mellor  et  al.  (1998),  the  present 
exclusion  of  the  prompt  mechanism  at  Diesel  conditions  is  further  supported  by  Dec 
(1997).  Using  laser-sheet  imaging,  he  shows  that  premixed  combustion  in  DI  Diesels  is 
not  conducive  to  NO  formation  through  the  prompt  mechanism  since  (|)  «  4.  Dec  (1997) 
does  argue  that  prompt  NO  may  form  in  the  diffusion  portion  of  the  Diesel  burning 
process;  nevertheless,  given  the  temperature  of  Diesel  combustion,  NO  formation  through 
this  pathway  is  likely  to  be  negligible  compared  to  the  alternatives. 

It  is  also  important  to  note  that  (R5)  is  a  three-body  reaction  and  is  thus  pressure 
sensitive.  In  the  following  section,  it  will  be  shown  that  this  pressure  sensitivity  changes 
the  relative  contribution  of  the  N2O  pathway  to  NO  formation  in  Diesel  engines. 

1.3  TWO-ZONE  MODEL  FOR  DIESEL  COMBUSTION 

Based  on  recent  laser  sheet  imaging  results  of  Dec  (1997)  and  the  earlier  findings 
of  Ahmad  and  Plee  (1983),  Mellor  et  al.  (1998)  have  postulated  a  two-zone  model  for  NO 
formation  and  decomposition  in  DI  Diesel  engines  as  illustrated  in  Fig.  1-1.  In  the  figure. 
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Figure  1-1.  The  two-zone  model  for  NO  formation  and  decomposition  in  a  DI  Diesel 
combustion  chamber  as  postulated  by  Mellor  et  al.  (1998).  Zone  1'  is  intermediate 
and  based  on  the  zone  1  pressure  and  flame  temperature  but  the  zone  2 
stoichiometry,  <|)ov!  for  the  engine. 
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Zone  1  is  a  near  stoichiometric  region  ((j) «  1)  where  [NO]  is  small  and  in  which  NO 
formation  occurs  through  forward  reactions  (Rl)  through  (R5),  equilibrium  (R6),  and 
forward  reaction  (R7).  Following  the  results  of  Dec  (1997),  this  zone  is  shown  as  a  thin 
layer  surrounding  the  periphery  of  the  fuel  plume.  Between  zone  1  and  zone  2,  there  is  a 
transition  in  which  the  concentrations  of  the  products  in  zone  1  are  adjusted  to  the  overall 
engine  stoichiometry  (zone  1').  Zone  2  is  a  region  in  which  hot,  burned  gases  are  present 
at  the  overall  equivalence  ratio.  In  this  zone,  NO  decomposition  occurs  through  the 
reverse  of  the  mechanism  given  above. 

The  adiabatic  stoichiometric  flame  temperature  T4,=i  for  zones  1  and  T  is 
computed  following  Plee  et  al.  (1980),  that  is,  with  start  of  combustion  (SOC)  conditions, 
the  stoichiometric  amount  of  fuel,  and  assuming  constant  pressure  burning. 
Stoichiometric  temperatures  vary  from  2200  K  to  2700  K  with  the  lower  values  resulting 
from  EGR  or  water  injection.  The  representative  temperature  in  zone  2  is  the  end  of 
combustion  temperature.  Tab  and  is  calculated  following  Mellor  et  al.  (1998).  Typical 
values  for  Tab  are  in  the  range  of  1300  to  2200  K  depending  upon  overall  equivalence 
ratio,  inlet  manifold  conditions,  and  so  forth. 

In  the  zonal  model  just  described,  it  is  assumed  that  concentration  and 
temperature  gradients  can  be  ignored  within  a  given  zone.  Thus,  Mellor  et  al.  (1998) 
model  the  chemistry  in  each  zone  using  a  perfectly  stirred  reactor  (PSR)  simulation 
(Bowman  et  al.,  1995)  with  adiabaticity  assumed.  The  GRI MECH  2.1 1  methane-air 
kinetic  mechanism  (Bowman  et  al.,  1995)  is  used  to  provide  heat  release  and,  in  this 
sense,  model  stoichiometric  and  lean  Diesel  fuel  combustion. 
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The  relatively  constant  concentrations  of  N-atom  and  N2O  in  each  zone  provide  a 
basis  for  the  steady-state  approximation  to  be  made  in  the  following  section.  Also,  the 
PSR  simulation  of  Mellor  et  al.  (1998)  shows  that  [NO]eq  is  largest  in  zone  1'  where  ^  = 
(j)ov  and  the  reactor  temperature  and  pressure  remain  at  the  stoichiometric  zone  conditions. 
This  fact  will  be  used  in  the  definition  of  the  characteristic  times  derived  below. 

1.4  NO, FORMATION  AND  DECOMPOSITION  TIMES 
1  A.  1  Characteristic  Time  For  NO  Formation 

Writing  the  law  of  mass  action  for  the  rate  of  NO  formation  in  zone  1  gives 
d[NO]/dt  =  kif[0][N2]  +  k2f[N][02]+  k3f[N][OH]  +  2k4f[N20][0]+  2k7f[N20][H]  (1) 

Invoking  the  N-atom  and  N2O  steady  state  assumptions  as  discussed  above,  Eq.  (1) 
becomes 

d[N0]/dt  =  2[0][N2](k,f+k5f)  (2) 

The  NO-forming  eddies,  to  which  Eq.  (2)  applies,  originate  just  inside  zone  1  and  are 
convected  outward  through  zone  1'  to  zone  2. 

Mellor  et  al.  (1998)  let  Sno  represent  the  fractional  extent  of  reaction  of  NO  such 
that  Sno  =  [NO]/[NO]eqr  where  0  <  sno  <  1  and  [NOjeqr  is  the  equilibrium  NO 
concentration  from  the  zone  1  stoichiometric  flame  temperature  calculation  adjusted  to 
the  overall  engine  stoichiometry.  As  shown  by  Mellor  et  al.  (1998),  [NOjeqr  is  the  proper 
choice  as  a  normalizing  factor  to  constrain  8no  between  the  limits  of  zero  and  one  since  it 
is  the  largest  value  among  the  various  zones. 

Substituting  this  definition  into  Eq.  (2)  gives 

dsNo/dt  =  (l/[NO]eqiO(ci[NO]/dt)  =  2[0][N2](kif  +  k5f)/[NO]eqi'  (3) 

The  inverse  of  this  equation  is  the  characteristic  time  for  NO  formation 
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t„o  =  [NO]eqr  /  (2[0]eql[N2]eql  (kjf  +  kjf ))  (4) 

where  [0]eqi,  [N2]eqi,  and  the  rate  coefficients  are  evaluated  at  the  zone  1  burned  gas 
conditions  (i.e.,  the  output  from  the  stoichiometric  equilibrium  calculation  based  on 
motored  engine  conditions  at  start  of  combustion). 

Also,  Mellor  et  al.  (1998)  define  Zeldovich  and  N2O  mechanism  formation  times 
as  follows: 

[NO]eql72klf[0]eql[N2]eql  (5) 

tno,5f=  [NO]eql72k5f{0]eql[N2]eql  (6) 

The  activation  temperatures  for  Tno,5f  and  Xno.ir  are  36120  K  and  62184  K,  respectively. 
The  latter  value  is  typical  of  the  values  obtained  for  NOx  emissions  from  conventional 
gas  turbine  combustors  (see  the  discussion  above). 

Recall  that  kif,  in  Eq.  (4)  and  (5),  represents  the  extended  Zeldovich  mechanism, 
(Rl)  through  (R3),  while  the  ksr  term  in  Eq.  (4)  and  (6)  is  the  third-body  component,  thus 
sensitive  to  pressure,  originating  in  the  N2O  mechanism,  (R4),  (R5),  and  (R7).  In  Fig.  1- 
2  each  contribution  in  Eq.  (2)  is  evaluated  in  Arrhenius  format  for  typical  engine 
conditions  using  the  kinetic  information  presented  in  Mellor  et  al.  (1998). 

The  figure  shows  that  at  low  pressures,  1  or  10  atm,  NO  formation  will  follow 
extended  Zeldovich  kinetics  (kif  exceeds  ksf)  over  most  of  the  temperature  range  shown. 
However,  at  10  atm  and  the  lower  stoichiometric  flame  temperatures,  obtained  with  inert 
injection  for  example,  the  N2O  mechanism  will  be  predominant  (ksf  exceeds  kif).  The 
crossover  temperature  moves  to  somewhat  higher  values  as  the  pressure  is  increased,  and 
for  60  atm  (representative  of  motored  TDC  Diesel  conditions),  the  N2O  mechanism 
dominates  at  many  stoichiometric  temperatures  of  interest  (approximately  2200  to  2460 
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Figure  1-2.  The  contribution  of  the  Zeldovich  and  NjO  mechanisms  to  the  NO 
formation  rate  as  postulated  by  Mellor  et  al.  (1998).  The  information  in  this  figure 
is  presented  in  Arrhenius  format  for  typical  engine  conditions.  Logarithms  are 
presented  as  functions  of  inverse  stoichiometric  flame  temperature  in  the 
approximate  range  from  2200  to  2900  K,  so  that  midrange  values  correspond  to 
stoichiometric  values,  with  pressure  as  a  parameter. 
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K).  The  contribution  of  this  mechanism  at  higher  values  of  stoichiometric  temperature 
can  not  be  ignored. 

1.4.2  Characteristic  Time  for  NO  Decomposition 

The  characteristic  time  for  NO  decomposition  is  defined  in  a  manner  analogous  to 
the  characteristic  time  for  NO  formation.  Writing  the  law  of  mass  action  for  the  rate  of 
NO  decomposition  in  zone  2  gives  (for  reasons  given  by  Mellor  et  al.  (1998),  the  reverse 
of  (R7)  is  ignored  for  the  present  analysis) 

d[NO]/dt  =  -  k,r[NO][N]  -  k2r[N0][0]  -  k3r[NO][H]  -  2k4r[NO]^  (7) 

Expressed  in  terms  of  the  fractional  extent  of  reaction,  and  again  assuming  steady-state 
N-atom,  this  equation  becomes 

dsNo/dt  =  (l/[NO]eqi')(d[NO]/dt)  =  (-  2kir[NO][N]  -  2k4r[NO]")/[NO]eqr  (8) 
The  inverse  of  this  expression  is  the  characteristic  time  for  NO  decomposition 

^no,3b  1/  [EnO  (2kii- [N]eq2  2k4|-[NO]eql'SNo)l  (9) 

The  form  of  this  equation  demonstrates  that  NO  decomposition  is  a  function  of  both  the 
reverse  Zeldovich  mechanism  (2kir[N]eq2)  and  the  reverse  N2O  mechanism 
(2k4r[NO]eqi'8No)-  Individual  decomposition  times  for  the  extended  Zeldovich  and  N2O 
mechanisms  are  respectively  expressed  as 


'tab.lr  -  l/2kir[N]eq2 

(10) 

’t3b,4r  l/2k4r[NO]eqr 

(11) 

In  an  Arrhenius  plot  similar  to  Fig.  1-2,  Mellor  et  al.  (1998)  show  that  NO  decomposition 
through  the  reverse  N2O  mechanism  is  dominant  over  a  range  of  typical  DI  Diesel 
operating  conditions  (p  =  60  atm  and  Tab  =  1300  K  to  2200  K). 
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1.5  CHARACTERISTIC  TIME  MODEL  FOR  NOx  FROM  DIESEL  ENGINES 
Using  the  kinetic  times  defined  above,  expressions  for  the  engine-out  NOx 
emissions  are  obtained  by  integrating  the  laws  of  mass  action  over  the  fluid  or  residence 
times  available  in  zones  l-T  and  2.  For  NO  formation  in  zones  1-1',  Eq.  (3)  and  (4)  are 
combined  such  that 

dSNo/dt  =  1/Tno  (11) 

where  [O]  and  [N2]  in  Xno  are  evaluated  at  constant  equilibrium  zone  1  conditions. 
Equation  (11)  can  be  integrated  over  the  eddy  lifetime  in  zone  1-1',  Tf,,no,  to  obtain 

SNor  =  ni'(Tr,.„o  /t„J  =  m’Da„„  (12) 

In  Eq.  (12),  8Nor  is  the  fractional  extent  of  reaction  for  NO  formed  in  zone  1-1'  per 
stoichiometric  eddy,  m'  is  a  scale  factor  introduced  because  the  magnitude  of  the 
characteristic  time  ratio  is  unknown  at  this  point,  and  Dano  is  the  Damkohler  number  for 
NO  formation  in  zone  1-1'. 

If  a  portion  of  the  NO  formed  in  zone  1-1'  decomposes  in  zone  2,  then  as  Mellor 
et  al.  (1998)  have  shown  in  detail,  Eq.  (8),  (10),  and  (11)  can  be  used  to  obtain  engine-out 
emissions  such  that 


NO,,EI  = 


m 


1  +  8 


Nor 


Da„„[NO],,,e-^’*- 


k4r[NO]^,,/ 


(13) 


where  m  is  a  model  constant  incorporating  the  scaling  factor  m'  and  other  conversion 
constants,  Dasbjr  is  a  Damkohler  number  based  on  the  eddy  lifetime  in  zone  2  (Tsi,3b)  and 
the  reverse  Zeldovich  mechanism  for  zone  2  (Dasb  =  Tsi,3b/x3b,ir),  nif  is  the  fuel  supplied 
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per  cycle  in  g,  and  Vab  is  the  volume  in  the  cylinder  at  the  end  of  combustion  m  cm  . 

Mellor  et  al.  (1998)  show  that  Eq.  (13)  can  be  simplified  for  various  limiting 
cases.  In  particular,  if  the  eddy  lifetime  in  zone  2  is  inadequate  so  that  NO 
decomposition  is  negligible  (i.e.,  Tspb  and  Daab.ir  approach  zero),  then 

NO,EI  =  m^Da„[NO]^,,  (14) 

mj. 

Another  limiting  case  of  interest  consists  of  adequate  time  for  zone  2  decomposition  with 
dominant  N2O  kinetics  (kir  ->  0),  the  most  likely  case  for  decomposition  as  shown  by 
Mellor  et  al.  (1998).  For  this  situation,  Eq.  (13)  reduces  to 

NO,EI  =  m^Da3t4.[NOL„,  (15) 

m, 

where  the  kinetic  time  in  Da3b,4r  is  evaluated  at  Tab. 

1.6  APPLICATION  OF  NOx  KINETIC  MODEL  TO  DI  DIESELS 
1.6.1  Determining  Significance  of  NO  Decomposition  in  Diesel  Engines 

The  ratio  of  the  characteristic  times  for  NO  decomposition  and  formation, 
Xab.no/'tno,  provides  a  method  to  evaluate  the  importance  of  the  former  in  any  engine 
situation.  This  ratio  is  developed  and  verified  in  Mellor  et  al.  (1998).  The  advantages  of 
using  this  ratio  are  its  ease  of  calculation  and  that  it  depends  only  on  the  fiiel  and  engine 
operating  parameters;  engine  load,  pressure  at  SOC  (P2),  temperature  at  SOC  (T2), 
compression  ratio  (rc),  and  maximum  cylinder  pressure  (Pa). 

The  effect  of  various  engine  operating  parameters  on  the  decomposition  to 
formation  time  ratio  is  investigated  by  conducting  an  equilibrium  analysis  of  various 
representative  engine  operating  conditions  for  four  DI  Diesel  engines.  The  engines 
examined  (three  heavy  duty  and  one  light  duty)  are  relatively  new  engines,  some  of 
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which  are  equipped  with  state-of-the-art  NOx  control  systems  (EGR  and  water  injection). 
Configurations  with  turbocharging  and  inter/aftercooling  as  well  as  naturally  aspirated 
are  available.  A  complete  table  of  engine  descriptions  is  given  in  Mellor  et  al.  (1998). 

The  variety  of  engine  setups  gives  a  wide  range  of  operating  conditions  from  which  to 
examine  the  results  of  the  study. 

Results  of  this  investigation  lead  to  the  following  conclusions  about  the  effects  of 
engine  operating  conditions  on  X3b,no/Xno-  By  far  the  most  important  engine  operating 
parameter  affecting  NO  decomposition  is  load.  The  higher  cylinder  temperatures  and 
pressures  associated  with  higher  loads  increase  the  rate  of  decomposition.  SOC  and  peak 
pressures,  injection  timing,  and  fuel  properties  have  negligible  effects.  Start-of- 
combustion  temperature  and  compression  ratio  have  small  influences  masked  by  changes 
in  load.  Inert  injection  results  can  be  interpreted  in  terms  of  their  effect  on  equilibrium 
concentrations  and  characteristic  temperatures  in  each  zone. 

1 .6.2  Fluid  Mechanic  Mixing  Times 

The  CTM  equation  for  the  prediction  of  engine-out  NOxEI  is  given  in  Eq.  (13).  In 
order  to  use  this  equation  values  for  the  fluid  mechanic  mixing  times  tf,,no  and  Tsi,3b  must 
be  known.  A  characteristic  mixing  time  can  be  defined  as  the  time  an  eddy  resides  in  a 
characteristic  zone  before  its  temperature  is  lowered  to  the  point  that  the  chemistry  of 
interest  is  quenched  by  dilution  or  expansion.  A  primary  assumption  of  the  model 

\ 

construction  done  thus  far  has  been  that  NO  forms  primarily  in  zone  one  and  decomposes 
in  zone  two.  Therefore,  the  characteristic  fluid  mechanic  mixing  time  for  NO  formation 
(Xfi,no)  is  defined  as  the  residence  time  of  a  stoichiometric  eddy  in  zone  1.  Likewise,  the 
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characteristic  mixing  time  for  NO  decomposition  (Xsi,3b)  is  defined  as  the  residence  time 
of  an  eddy  in  zone  2. 

Equation  (13)  includes  both  formation  and  decomposition  reactions.  Of  the 
three  unknowns  in  the  equation,  two  are  mixing  times  and  one  is  a  model  constant.  The 
model  constant  is  equal  to  the  slope  of  the  linear  best-fit  line  formed  by  graphing  engine- 
out  NOxEI  versus  the  appropriate  model  parameter.  Thus,  the  constant  can  be  solved  for 
once  the  other  two  unknowns  have  been  determined  and  is  therefore  ignored  in  the 
following  analysis.  The  characteristic  mixing  time  for  NO  decomposition  is  eliminated 
fi'om  the  equation  by  utilization  of  low  load  engine  data  for  which  NO  decomposition  is 
thought  to  be  negligible  (Tno,3b»  Xno),  thereby  simplifying  the  equation  to  Eq.  (14)  with 
only  Tf,,no  unknown. 

Engine  C  data  are  used  for  the  tfi,no  study  because  these  data  contain  a  wide  range 
of  operating  conditions  at  low  load  ('Cno,3b»Xno).  Engine  C  is  a  HSDI  Diesel  with  a 
displacement  of  0.55  L/Cyl  and  is  equipped  with  a  turbocharger  and  intercooler.  Mellor 
et  al.  (1998)  give  a  full  description  of  the  engine.  Analysis  of  Engine  C  data  indicates 
that  NO  decomposition  can  be  assumed  negligible  for  loads  of  600  kPa  and  below  as 
discussed  in  Mellor  et  al.  (1998). 

An  ongoing  literature  review  gives  insight  as  to  which  engine  parameters  should 
be  included  in  the  mixing  time.  Some  of  the  topics  being  studied  along  with  key 
references  are  listed  in  Table  1-1 .  This  list  is  by  no  means  thorough,  but  does  give  an 
idea  of  the  mixing  processes  examined  in  the  study. 

Knowledge  gained  from  the  literature  review  about  in-cylinder  mixing  is  used  to 
compile  a  list  of  mixing  processes  important  to  fuel-air  mixing  in  the  cylinder  and  the 
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characteristic  terms  by  which  they  are  modeled.  Some  such  terms  are:  Weair  -  secondary 
breakup,  Reair  -  turbulence  in  air,  etc.  Using  this  list  as  a  starting  point,  different 


Table  1-1.  Fluid  Mechanic  Research  Topics 


Effects  of  cylinder  gas  density  on  spray 
penetration  and  dispersion 

Nabers  and  Siebers  (1996),  Arai  and 
Hiroyasu  (1990) 

Fuel  film  formation  and  wall  interaction 

Stanton  and  Rutland  (1996) 

Engine  speed 

Uludogan  et  al.  (1996) 

Sauter  mean  diameter 

Hiroyasu  et  al.  (1989) 

Secondary  breakup 

Arcoumanis  et  al.  (1997) 

Effects  of  swirl 

Timoney  and  Smith  (1996) 

Overall  combustion  and/or  mixing  sources 

Arcoumanis  et  al.  (1997),  Heywood 
(1988),  Ferguson  (1986) 

definitions  of  the  mixing  time  are  tested  to  establish  the  best  correlation  between  engine- 
out  NOxEI  and  the  CTM  model  parameter. 

A  single  tfi,no  for  all  Engine  C  data,  assuming  NO  decomposition  negligible,  has 
not  yet  been  obtained.  Instead,  the  data  are  correlated  for  the  present  using  two  separate 
empirical  mixing  times  for  NO  formation.  The  first  applies  to  cases  in  which  fuel 
impinges  on  the  cylinder  bowl  before  SOC  (wall  jet),  and  the  second  to  cases  in  which 
fuel  impinges  after  SOC  or  not  at  all  (free  jet).  The  mixing  times  contain  identical  terms 
that  are  related  to  spray  breakup,  impingement  time,  and  ignition  delay  time,  but  the 
significance  (i.e.  exponent)  of  each  term  is  different  for  the  two  cases.  The  mixing  time 
for  the  wall  jet  case  is  defined  as 


We'  ®  1  d  T 

_  *imp  '^noz  ✓ ''imp  v -0.5 

^fi.no  5  i  I 

'“'fuel  ‘’id 


and  the  mixing  time  for  the  free  jet  case  is  found  to  be 


(16) 
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Graphs  of  NOxEI  versus  the  right  hand  side  of  Eq.  (14)  can  be  found  in  Fig.  1-3 
and  1-4.  These  graphs  use  data  that  meet  the  criteria  of  negligible  decomposition  (load  < 
600kPa)  and  SOC  =  0.  The  linear  best-fit  lines  give  correlation  coefficient,  R,  values  of 
0.9793  and  0.9704  for  the  impingement  before  SOC  and  impingement  after  SOC  cases, 
respectively.  A  problem  arises  with  these  linear  fits  in  that  as  Xno  goes  to  infinity  the 
engine-out  NOx  is  expected  to  go  to  zero,  but  for  both  cases  the  y-intercept  is  relatively 
large.  We  are  evaluating  the  method  of  defining  Xno  at  this  time.  The  lowest  engine-out 
NOx  datum  points  seem  to  indicate  the  expected  y-intercept  of  zero,  however.  The 
graphs  presented  are  preliminary  and  appropriate  fits  will  be  determined  after  further 
research. 

Currently,  our  efforts  are  directed  toward  improving  the  empirical  mixing  times 
for  NO  formation.  We  are  hopeful  that  a  single  mixing  time  for  NO  formation  can  be 
found  that  is  valid  for  start  of  combustion  prior  to  and  following  wall  impingement. 

Once  a  correctly  defined  mixing  time  for  NO  formation  is  developed  for  engine  C,  it  will 
be  validated  for  other  engines.  With  a  satisfactory  definition  of  the  mixing  time  for  NO 
formation,  work  will  begin  on  the  mixing  time  for  NO  decomposition.  When  both 
mixing  times  have  been  defined  the  CTM  will  be  capable  of  predicting  NOxEI  for  any  DI 
Diesel  operating  condition. 

1.7  CHARACTERISTIC  TIME  MODEL  FOR  SOOT  FROM  DIESEL  ENGINES 

Although  the  model  for  soot  emissions  is  in  the  early  stages  of  development,  it  is 
postulated  that  a  two-zone  model  representing  the  formation  and  oxidation  of  soot  will 
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Figure  1-3.  Engine-out  NOx  emissions  index  versus  model  parameter  (Eq.  (14))  for 
wall  jet  datum  points  with  SOC  =  0.  The  legend  has  the  format:  engine  speed 
(RPM)/load  (kPa)/injection  pressure  (MPa).  Each  set  of  data  is  obtained  via  an 
EGR  sweep. 
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Figure  1-4.  Engine-out  NOx  emissions  index  versus  model  parameter  (£q.  (14))  for 
free  jet  datum  points  with  SOC  =  0.  The  legend  has  the  format:  engine  speed 
(RPM)/load  (kFa)/injection  pressure  (MPa).  Each  set  of  data  is  obtained  via  an 
EGR  sweep. 
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properly  model  this  exhaust  emission.  Figure  1-5,  based  on  past  laser  sheet  images  (see 
Dec,  1997),  shows  a  schematic  of  a  Diesel  spray  plume  with  the  pertinent  zones  where 
soot  is  forming  and  oxidizing.  Dec  (1997)  postulates  that  during  the  diffusion  flame 
(mixing  controlled)  combustion,  a  thin  premixed  flame  will  form  at  the  liquid/vapor  fuel 
interface,  and  initiation  of  soot  occurs  just  downstream  of  this  flame  in  the  products  of 
rich  combustion.  Soot  will  then  continue  to  form  in  the  vapor  fuel  core,  where  a 
hypothetical  homogeneous  mixture  with  (|)»3  exists,  until  it  reaches  the  diffusion  flame 
front.  Thus,  soot  formation  is  assumed  to  occur  spatially  in  the  vapor  core  of  the  spray, 
labeled  zone  1  in  Fig.  1-5.  For  the  oxidation  of  soot,  Dec  (1997)  shows  that  due  to  the 
large  concentrations  of  hydroxyl  (OH)  radicals  in  the  diffusion  flame,  the  majority  of  soot 
oxidation  will  occur  in  this  location.  This  assumption  is  further  supported  through  the 
work  of  Dec  and  Coy  (1996)  and  Kosaka  et  al.  (1995).  For  modeling  purposes,  it  will  be 
assumed  that  soot  oxidation  is  initiated  at  the  stoichiometric  contour  of  the  fuel  spray  and 
continues  out  into  the  periphery  of  the  flame  (see  Fig.  1-5). 

The  representative  temperature  of  zone  1 ,  Tsf,  is  defined  as  an  average  of  2 
limiting  temperatures,  namely  the  50%  fuel  recovery  temperature,  T5o%,  and  the  (t)=3 
flame  temperature,  T,|,=3,  where  T,j,=3  is  calculated  using  TDC  cylinder  conditions  as  done 
for  T,|^i.  The  temperature  indicative  of  zone  2,  Tso,  is  defined  as  the  average  of  the 
stoichiometric  and  end  of  combustion  temperatures,  T^=i  and  Tsb. 

Assuming  the  presence  of  soot  formation  and  oxidation,  the  emissions  model  for 
soot  was  developed  following  the  approach  used  by  Mellor  et  al.  (1998)  to  derive  the 
NOx  CTM,  except  that  global  reactions  for  soot  formation  and  oxidation  will  be  used 
instead  of  a  skeletal  mechanism  because  of  the  complex  chemistry  involved  in  the 
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Figure  1-5.  Schematic  of  a  DI  Diesel  spray  plume  with  regions  of  soot  formation 
and  oxidation  labeled. 
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formation  and  oxidation  of  soot.  This  is  consistent  with  the  work  of  others  to  model  soot 
emissions  from  Diesel  engines  (see  for  example,  Nishida  and  Hiroyasu,  1989). 

Following  the  methodology  used  in  developing  the  NOx  CTM,  the  overall 
reactions  used  to  model  soot  formation  and  oxidation  are  as  follows; 

Fuel —>  Soot  (18) 

Soot  +  OH^C02  (19) 

Applying  the  law  of  mass  action  to  Eqs.  (18)  and  (19)  gives 

=k^[fuel]  (20) 

dt  form 

=  -k,Jsoot].[OH]  (21) 

oxid 

In  Eq.  (20)  it  is  assumed  that  soot  has  a  first  order  of  dependence  on  fuel.  Now 
integrating  Eqs.  (20)  and  (21)  over  the  lifetime  that  eddies  spend  in  the  soot  formation 
and  oxidation  zones,  Tsi,sf  and  Xsi,so,  respectively,  yields 


d[soot] 

dt 


[soot],  _ 
[fuel] 


=  m^Da^ 


where  Tsf=  1/ksf,  Tso=  l/kso[OH],  and  msf  and  mso  are  model  constants.  For  Eq.  (23)  it  is 
assumed  that  [OH]  is  in  steady  state,  which  is  based  on  the  work  of  Tuttle  et  al.  (1977) 
for  the  CO  CTM.  However,  further  development  of  this  model  by  Connors  et  al.  (1996) 
incorporates  a  first  order  dependence  of  [OH].  The  concentrations  [soot]i  and  [soot]2  in 
Eq.  (23)  represent  the  amount  of  soot  formed  in  zone  1  and  the  total  soot  remaining  after 
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oxidation  in  zone  2  has  ceased,  where  the  amount  [soot]2  is  dependent  on  the  soot  formed 
in  zone  1 . 

The  kinetic  times  associated  with  each  of  the  chemical  processes  will  be 
represented  as  inverse  Arrhenius  expressions: 

X,  -  exp(E  JRTJ  (24) 

with  Ex  and  Tx  the  global  activation  energy  and  representative  temperature,  respectively, 
for  each  process. 

The  engine-out  soot,  in  units  of  emissions  index  of  soot,  is  obtained  by  combining 
Eqs.  (22)  and  (23): 

SootEI  oc  =  m,f Da,f  •  e"”"''*-”  (25) 

[fuel] 

The  soot  model,  Eq.  (25),  includes  the  competition  between  soot  formation  and  oxidation 
with  each  process  represented  by  the  Damkohler  numbers  Dasf  and  Dso,  respectively. 

Tuttle  et  al.  (1977)  showed  the  exponential  term  in  the  CO  CTM,  similar  to  that  in 
Eq.  (25),  can  be  approximated  by  a  series  expansion  of  Tso/Tsi,so: 

+  ...  (26) 

where  the  coefficients  ao,  ai,  a2  . . .  are  determined  through  best  fitting.  Furthermore,  Tuttle 
et  al.  (1977)  found  that  neglecting  the  squared  and  higher  terms  for  their  measurements 
results  in  errors  of  less  than  one  percent.  If  the  same  holds  true  for  soot  data,  then  the 
exponential  term  in  Eq.  (25)  can  be  represented  by  the  two  term  expansion  resulting  in  the 
following: 


SootEI  oc  m^Da^f 


(27) 
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Equation  (27)  represents  the  final  soot  model,  where  the  coefficient  nisf  is  determined  from 
the  best-fit  slope  of  experimental  data  of  SootEI  plotted  versus  Dasf  times  the  bracketed 
term.  Although  Eq.  (27)  is  preliminary,  a  point  of  departure  would  be  to  assure  that  the 
combination  of  both  the  NOx  and  soot  models  is  capable  of  predicting  the  trends  existing 
between  these  exhaust  emissions  under  various  engine  operating  conditions. 

1.8  NOx  -  SOOT  EMISSION  TRENDS 
1.8.1  Diesel  Engine  Performance 

Since  the  Diesel  engine  is  used  in  a  number  of  applications,  fi-om  large,  heavy- 
duty  (HD)  engines  for  the  military  in  combat  vehicles  to  the  smaller  light-duty  (LD) 
Diesels  in  small  trucks  and  automobiles,  performance  requirements  demand  engine 
operation  over  wide  load  and  speed  cycles.  The  transient  conditions  under  which 
commercial  vehicles  operate  require  Diesel  engines  to  typically  run  at  three  quarter  to  full 
load  (see  Mori,  1997),  while  military  vehicles  spend  the  majority  of  time  at  low  load 
conditions.  These  different  load  cycles  result  in  different  engine  performance 
requirements  which  in  turn  lead  to  changes  in  the  trends  existing  between  NOx  and  soot 
emissions. 

Since  any  emissions  model  must  be  capable  of  predicting  the  changes  in  trends  of 
exhaust  emissions  with  varying  engine  operating  conditions,  examination  of  a 
performance  map  for  DI  Diesel  engines  is  beneficial.  Figure  1-6,  atypical  performance 
map,  describes  the  effects  of  load  and  speed  variations  on  bsfc  (Heywood,  1988).  The 
performance  map  shows  that  at  approximately  50%  load  and  speed,  the  best  fuel 
economy  of  the  engine  is  obtained.  Note  for  the  engine  of  Fig.  1-6,  the  maximum 
attainable  load  is  limited  to  about  800  kPa  bmep  due  to  the  smoke  limit  of  the  engine. 
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bmep,  kPa 


Mean  piston  speed,  m/s 


4  6  8  10  12 


Figure  1-6.  Performance  map  for  a  naturally  aspirated  DI  Diesel  engine  (Heywood, 
1988). 
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Torque,  N*m 


The  Blue  Ribbon  Committee  Report  (BRC,  1995)  shows  the  tradeoff  between 
load  (bmep)  and  bsfc  through  the  relationship  to  the  power  density  of  a  Diesel  engine 
(pd),  defined  as  the  engine  output  divided  by  the  volume  of  the  engine  plus  accessories 
(i.e.,  turbocharger,  air  filter,  etc.): 

Pd  oc  bmep  oc(pa)((j))(l/bsfc)  (28) 

where  pa  is  the  cylinder  inlet  air  density  and  (j)  the  equivalence  ratio.  From  Eq.  (28),  the 
highest  Pd  would  be  obtained  by  a  turbocharged  (increased  pa)  engine  operating  under 
rich  conditions  with  low  bsfc;  however,  rich  operating  conditions  will  result  in  higher 
particulate  and  bsfc  (BRC,  1995).  The  tradeoff  in  Eq.  (28)  between  bmep  and  bsfc 
translates  into  different  trends  existing  between  NOx  and  soot. 

The  work  of  Wilson  et  al.  (1973)  and  Murayama  et  al.  (1978)  clearly  exhibits  the 
effects  that  increasing  load  has  on  bsfc  and  emissions  of  soot  and  NOx  for  a  constant 
speed.  As  the  load  (or  (|))  is  increased,  NOx  and  soot  emissions  simultaneously  increase 
until  a  maximum  value  of  NOx  is  reached,  after  which  a  further  increase  in  load  yields  a 
decrease  in  NOx  and  a  continued  increase  in  soot.  The  initial  increase  in  NOx  with  load 
corresponds  to  the  decrease  in  bsfc;  once  the  minimum  in  bsfc  is  obtained,  roughly 
coinciding  with  the  maximum  NOx,  any  further  increase  in  load  results  in  a  decrease  in 
NOx  that  is  caused  not  only  by  the  increase  in  bsfc,  but  also  extended  burning  into  the 
expansion  stroke  causing  reductions  in  peak  gas  temperature  and  increased  heat  losses  so 
the  production  rate  of  NO  is  reduced,  according  to  Wilson  et  al.  (1973).  As  discussed  in 
the  previous  sections,  an  alternative  interpretation  of  the  results  involves  increased 
decomposition  of  NO  as  load  is  increased. 
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The  above  discussion  suggests  that  there  are  two  different  trends  present  between 
NOx  and  soot  emissions  with  varying  load:  1)  low  load  cases  (<  50%  load)  show  an 
increase  in  NOx  and  soot  with  increasing  load  and  2)  higher  load  cases  (>50%)  where  the 
typical  NOx/soot  tradeoff  exists.  Algegraic  expressions  are  derived  below  from  the  NOx 
and  soot  CTMs  to  prove  that  the  overall  model  will  properly  predict  the  trends  of  NOx 
and  soot  emissions,  regardless  of  operating  conditions. 

1.8.2  Derivation  of  NO^Soot  Trends  Using  the  CTM 

To  derive  the  NOx/soot  emissions  equations  for  Diesel  engines,  a  relationship 
between  the  NOx  and  soot  CTM  is  needed.  The  link  in  the  derivation  is  the  mixing  times 
involved  in  the  models  for  NOx  and  soot  emissions.  The  zonal  assumptions  of  the  NOx 
and  soot  CTMs  indicate  that  NO  is  forming  at  the  diffusion  flame  front,  which  as  noted 
previously  is  the  same  spatial  region  for  soot  oxidation  initiation.  Therefore,  a  possible 
relationship  between  these  two  mixing  times  may  exist  and  could  be  used  in  the 
derivation  of  the  NOx/soot  tradeoff  equations. 

Analysis  of  the  Engine  C  data^  leads  to  an  expression  for  the  mixing  time 
associated  with  NO  formation,  Tfi,no,  where  NOx  has  a  first  order  dependence  with 
pressure  drop,  AP,  across  the  injection  nozzle,  and  a  -1.6  order  of  dependence  with 
engine  speed  (Eq.  17).  A  similar  analysis  for  Engine  C  soot  data  was  performed  to 
determine  the  dependency  of  soot  with  each  of  the  above  operating  parameters;  however, 
the  chemical  kinetics  for  the  soot  model  are  not  yet  fully  developed,  so  engine  conditions 
must  be  such  that  only  the  engine  parameter  of  interest  (i.e.,  AP  or  N)  is  varying.  For 

^  Analysis  only  includes  data  where  NO  decomposition  is  negligible  (200-600  bmep  range)  and  the  injected 
fuel  ignites  prior  to  impinging  on  the  piston  bowl  wall  (free  jet). 
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example,  soot  data  for  constant  load,  EGR  and  injection  pressure  were  used  to  determine 
the  order  of  dependence  of  soot  with  engine  speed.  The  results  show  that  soot  has  an 
approximate  third  order  dependence  with  engine  speed,  which  is  twice  the  inverse  of  that 
found  for  NOx  and  slightly  higher  than  the  value  used  by  Plee  et  al.  (1981).  There  is  no 
single  order  dependence  of  soot  on  AP;  instead  the  magnitude  of  dependence  varies  with 
EGR  (i.e.,  increasing  EGR  decreases  the  order  dependence).  However,  in  general  the 
analysis  did  reveal  that  soot  emissions  will  have  an  inverse  dependency  on  AP,  opposite 
the  dependency  found  for  NOx,  which  is  expected. 

Several  major  issues  arose  during  the  above  analysis  of  the  soot  emissions  from 
Engine  C.  The  first  problem  was  that  the  yet-to-be-determined  chemical  kinetics  for  soot 
formation  and  oxidation  act  as  a  constraint  and  eliminate  many  datum  points  since  only 
AP  orN  variations  can  be  considered,  while  load,  EGR,  and  SOC  are  constant.  Thus, 
more  datum  sets  with  only  AP  or  N  changing  are  needed  in  order  to  ensure  that  the 
obtained  dependencies  of  soot  with  AP  and  N  are  well  understood.  Another  issue 
concerning  the  available  data  is  the  lack  of  high  load  run  conditions.  These  higher  load 
data  are  required  to  test  and  ensure  that  the  CTM  will  properly  predict  soot  over  the  entire 
load  cycle,  not  just  at  low  loads. 

Although  there  exist  some  questions  about  the  data  available  for  analysis,  the 
overall  trends  of  soot  emissions  with  AP  and  speed  are  consistent  with  those  reported  by 
Plee  et  al.  (1981);  therefore,  the  NOx/soot  proportionality  can  be  derived  for  these  data. 

In  general,  it  has  been  shown  that  soot  has  an  inverse  and  perhaps  a  higher  order  of 
dependence  on  AP  and  N  than  found  for  NOx.  These  results  would  suggest  the  mixing 
times  associated  with  soot  oxidation  and  NO  formation  are  inversely  proportional: 
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By  solving  for  Xfij,o  and  Tsi,so  in  Eqs.  (14)  and  (27),  respectively,  and  substituting  these  new 
expressions  into  Eq.  (29),  the  following  NOx/soot  equation  is  obtained: 


Equation  (30)  shows  that  the  NOx  and  soot  emissions  are  directly  proportional  for  the  lower 
load  data,  which  is  consistent  with  the  trends  shown  by  Wilson  et  al.  (1973)  and  Murayama 
et  al.  (1978).  Higher  load  data  will  be  required  to  determine  the  relationship  existing 
between  the  mixing  times  associated  with  NOx  and  soot  emissions.  The  higher  load  data  are 
expected  to  exhibit  the  NOx/soot  tradeoff  obtained  through  “normal”  vehicle  operation. 
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